Velocity Distribution of Topological Defects in Phase-Ordering Systems 
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The distribution of interface (domain-wall) velocities v in a phase-ordering system is considered. 
Heuristic scaling arguments based on the disappearance of small domains lead to a power-law tail, 
Pv(v) ~ v~ p for large v, in the distribution of v = |v|. The exponent p is given by p = 2 + d/ (z — 1), 
where d is the space dimension and 1/z is the growth exponent, i.e. z — 2 for nonconserved (model 
A) dynamics and z — 3 for the conserved case (model B). The nonconserved result is exemplified by 
an approximate calculation of the full distribution using a gaussian closure scheme. The heuristic 
arguments are readily generalized to systems described by a vector order parameter. 
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I. INTRODUCTION 



The theory of phase-ordering dynamics has seen signifi- 
cant advances in recent years pj] . Coarsening proceeds by 
the elimination of topological defects, e.g. domain walls, 
vortices, or strings, and various properties of the defect 
distribution have been investigated in some detail |^|,[| . 

In a recent paper, Mazenko Q| has carried out the first 
investigation of the distribution of defect velocities in a 
phase-ordering system. Using an approximate 'gaussian 
closure' scheme, he has computed the velocity distribu- 
tion for point defects in 0(n) models with nonconserved 
order parameter. He finds the interesting result that the 
velocity distribution has a power-law tail. 

In the present work we shall show that such power-law 
tails can be deduced rather generally using very simple 
scaling arguments. The central idea is that large defects 
velocities are associated with the vanishing of defects at 
small length scales. Hence the tail gives information 
about the small-scale structure of the defect morphol- 
ogy. However, the velocity distribution is more cleanly 
defined and therefore more convenient to study theoret- 
ically (and in simulations) than the morphology itself. 
The physical arguments we shall use are very general, 
and can be applied to systems with both point-like and 
extended defects, with or without conservation of the or- 
der parameter. 

Consider, for example, systems described by a scalar 
order parameter. This class includes the 'standard' 
phase-ordering systems such as binary alloys (order- 
disorder transitions or phase separation) and binary liq- 
uids. The configuration of the coarsening system is de- 
scribed by the locations of the interfaces, or domain walls, 
that separate regions occupied by the two ordered phases. 
We are interested in the velocity distribution of points on 
the interfaces. The high-velocity tail in the distribution 
arises as follows. At all stages of the coarsening process, 
small domains are being eliminated. When the domains 
are very small, the walls move very quickly. So the high- 
velocity tail is related to the density of small domains. 
This is the basic idea behind the scaling argument pre- 



sented in section II. For the case of non-conserved scalar 
fields, an explicit calculation using a 'gaussian closure' 
scheme (section III) gives an approximate form for the 
full velocity distribution, extending Mazenko's work Q] 
on the velocity distribution of point defects in 0(n) mod- 
els. The tail of the distribution agrees with the result 
derived from scaling arguments. 

We shall show that an alternative (and a priori equally 
valid) way of using the gaussian closure results to com- 
pute the velocity distribution gives a gaussian tail rather 
than a power-law tail. The defect velocity distribution, 
therefore, is a useful discriminator between various com- 
putational schemes. It is reassuring that, in the limit 
of large spatial dimension d, when the gaussian closure 
approximation is believed to become exact (for noncon- 
served fields), both schemes yield the same results. 

In section IV the scaling arguments are extended to 
vector fields, and the general result p — 2 + (d + I — 
n)/(z — 1) is obtained for the tail exponent. For point 
defects (n = d) this result agrees with that obtained by 
Mazenko Q using the gaussian closure approximation. 
On the basis of the results for scalar fields and vector 
fields with point defects, we conjecture [equation a 
general form for the defect velocity distribution within 
the gaussian approximation. 



II. SCALING APPROACH 
A. Nonconserved Scalar Fields 

A system described by a nonconserved scalar field 
coarsens by curvature driven growth [jlj. The normal 
velocity of a point on an interface is proportional to the 
total curvature K . For a small circular (or spherical, in 
d = 3) domain of radius r, K — (d — l)/r. The velocity 
of the interface is therefore v — dr/dt oc — 1/r, It fol- 
lows that the time t for a domain of initial radius r to 
disappear scales as t oc r 2 . 

Let n(r)dr be the number of domains per unit volume 
with radius between r and r + dr. The assumed scaling 
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of the domain morphology gives 

n{r) = L"( d + X ) f( r /L) 



(1) 



where L{t) ~ t 1 / 2 is the characteristic length scale at 
time t after the quench into the ordered phase. Integrat- 
ing (Q) over r gives of order one domain per scale volume 
L d , as required by scaling. 

In a time interval At, domains with r < (At) 1 / 2 will 
disappear. From (fil) the number of such domains is 



L -(d+l) 



(At) 1 



dr f{r/L) ~ L- (d+ V [dL/dt)M . (2) 



The right-hand side of follows from the fact that the 
total number of domains per unit volume scales as L~ d . 
The requirement that the left-hand side of @) be linear 
in At for At — > implies the property f(x) ~ x, for 
x —> 0, for the scaling function /(a;). Inserting this form 
in (J2J) , one recovers L(t) ~ t 1 / 2 as expected. 

We have shown that n(r)dr ~ rdr/L d+2 for r <C L. 
The interfacial area per unit volume associated with these 
small domains is 



^(r)dr - r d - 1 n{r) dr ~ r d dr/L d+2 



(3) 



Normalizing by the total interfacial area per unit volume, 
L , gives the area- weighted probability for interfacial 
radius of curvature between r and r + dr, i.e. 



P r (r)dr ~ r d dr / L d+1 , r < L 



(4) 



Since the velocity associated with radius-of-curvature r is 
v ~ 1/r, the interfacial velocity distribution P v (v) is ob- 
tained from P v (v) — P r {r)\dr / dv\, giving the power-law 
tail 



Pv(v) 



1 



v(vL) 



d+l 



(5) 



This, together with its generalizations below, is our main 
result. 



B. Conserved Scalar Fields 

For conserved scalar fields, coarsening pro- 
ceeds through the Lifshitz-Slyozov- Wagner (LSW) 
evaporation-condensation mechanism S , by which large 
domains (or, more generally, regions with low interfacial 
curvature) grow at the expense of small domains (regions 
of high curvature). This leads to the LSW growth law 
L(t) ~ t 1 / 3 . Large domain-wall velocities are associated, 
as in the non-conserved case, with the disappearance of 
small domains. For the conserved case, however, the 
relation between velocity and radius is v(r) ~ 1/r 2 
Therefore, the time taken for a domain of size r <C L(t) 
to evaporate is of order r 3 . In a time interval At, do- 
mains of size r < (At) 1 / 3 will disappear. Using ([!]), the 
number of such domains is 



L 



-(d+X) 



(At) 



1/3 



dr f(r/L) ~ L-( d+1 > {dL/dt) At , (6) 



in complete analogy to (0). The requirement that the 
left-hand-side be proportional to At forces the domain- 
size distribution function f{x) to have the small-x form 
f(x) ~ x 2 for conserved scalar fields. Inserting this form 
into @ one recovers L(t) ~ t 1 / 3 as expected. 

The interfacial area per unit volume associated with 



these small domains is given by A(r)dr 



„d-l 



n(r)dr, 



as in (||). Using n(r)dr ~ r 2 dr/L d+3 for r <C L, which 
follows from f(x) ~ x 2 , and normalizing by the interfa- 
cial area per unit volume, L _1 , gives the area- weighted 
probability for interfacial radius of curvature between r 
and r + dr: 



P r {r)dr - r d+i dr/L 



-d+2 



r < L 



(7) 



instead of (^). The final step is to use the relation 
v ~ 1/r 2 to deduce the velocity distribution P v (v) from 
P v (v) = P r (r)\dr / dv\. From <g) one infers the power-law 
tail 



Pv(v) 



v{vL 2 Y d + 2 )/ 2 



(8) 



for conserved scalar fields. 

The scaling approach adopted for scalar fields can be 
readily generalized to vector fields. Before doing so, how- 
ever, we will consider an analytical approach, based on 
a gaussian closure scheme, for computing the full distri- 
bution P v (v) (i.e. not just the tail) for the nonconserved 
scalar case. This is instructive as it provides, in con- 
junction with the scaling result @ for the tail, a way of 
discriminating between two (a priori equally valid) ways 
of using the gaussian closure results in the calculation of 
P, 



III. CALCULATION OF P V (V) USING GAUSSIAN 
CLOSURE 



The idea underlying all gaussian closure schemes is 
the introduction of a 'smooth' auxiliary field m(x, t) 
whose zeros define the interfaces, i.e. m has its zeros at 
the the same points as the order-parameter field 4>(x, t). 
Whereas tf> is essentially constant within domains, but 
varies rapidly within domain walls, m is 'smooth', i.e. it 
varies only on the larger scale Lit). At this stage, how- 
ever, it is not necessary to specify precisely how m is 
defined. 

An expression for the interface velocity in terms of 
m can be obtained by noting that the rate of change 
of m in a frame moving with the interface is zero, i.e. 
dm/ dt = = dm/dt + v • Vm. Since Vm is normal to 
the interface, v • Vm = w„|Vm|, where v n is the normal 
velocity of the interface. This gives 
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-dtml\ Vm| 



(9) 



Formally this equation defines a 'velocity' at every point 
- the velocity of the surface of constant m defined at that 
point. To find the distribution of interface velocities, we 
have to project onto the interfaces: 



Pv(v) 



(6(v - v n (x))p(x)) 
(P(x)) 



(10) 



where p(x) is the areal density of interface, given by 
p(x) = (5(m(x)) |Vm| . 



(11) 



Integrating p(x) over any volume of space gives the in- 
terfacial area in that volume. 

Equations - 11) are exact, but in order to make 
further progress one needs to know the joint distribution 
function for m, dtm, and the components d^m of Vm, at a 
given point in space. In all gaussian closure schemes, the 
field m(x, t) is assumed to be gaussian. This approxima- 
tion works reasonably well in practice for nonconserved 
fields [pj. Using it, the required distribution function is 
expressible in terms of the inverse of the covariance ma- 
trix of the variables m,m= dtm, and dim (i = 1, . . . , d). 
Clearly (mc^m) and (mc^m) both vanish due to transla- 
tional invariance. Similarly, (dimdjm) vanishes for i =/= j 
due to the assumed isotropy of the system. It follows 
that the only non-zero elements of the covariance ma- 
trix are S(t) = (m 2 ), 5/2 = (mm), T(t) = (m 2 ) and 
U(t) = ((dim) 2 ). Note that the Vm sector is decoupled 
from the (m, m) sector. The required probability distri- 
bution has, therefore, the product form 



P(m, m, Vm) 



(2wU) d / [ 



■ exp 



(Vm)< 
2U 



x ^rexp ( --^-(Tm 2 - 5mm + 5m 2 ) I , (12) 



2ttD 



where 



2D' 



D = ST- 574 



(13) 



is the determinant of the covariance submatrix in the 
(m, m) sector. 
From ( |Io| ) - ([ll]) the velocity distribution is give by 

o M _ ( S ( v + rh/\Vm\) S(m) |V(m)|) 
v{ ' (<5(m)|V(m)|) • 1 ' 

Evaluating the averages using the distribution ( |l2| ) gives 

1 rf^±2) 1 1 

v[ ) r(^) v (i + v yv 2 )( d + 2 )/ 2 ' [ ' 

where the characteristic velocity v is given by 




m 2 ) ((Vm) 2 )/d 



(16) 



Note that the distribution (15) is normalized on the in- 
terval (—00, 00). 

The y~( d + 2 ) tail predicted by ( ^5[ ) agrees with the scal- 
ing arguments for non-conserved fields presented in the 
previous section. This agreement is a nontrivial result, 
since up to now we have made no assumptions about 
the nature of the dynamics! The only assumption used 
is that the field m is gaussian, which is evidently quali- 
tatively correct for nonconserved fields but not for con- 
served fields, where the power-law predicted using (|l5| ) 
disagrees with the scaling result (||). We conclude that 
no simple gaussian auxiliary-field approach captures the 
physics of the conserved model. Similar conclusions have 
been drawn in another context || . 

The explicit value of v in ( |l6| ) depends on the dynam- 
ics of the field m. The simplest approximation is that of 
Ohta, Jasnow and Kawasaki (O JK) H , in which m obeys 
the diffusion equation, dtm — DoV^m,, with Do = (d — 
l)/d. The initial condition is (m(x, 0)m(x', 0)) = A<5(x— 
x'). With this choice one obtains v = [(d — l)/2i\ 1 / 2 . 
Other approximations, such as that of Mazenko M, give 
the same i -1 / 2 behaviour, but with a different prefactor. 
The £ -1 / 2 dependence follows quite generally from the 
scaling relation v ~ dL/dt and the usual result L ~ i 1 / 2 
for nonconserved scalar fields. 

Before concluding this section, it is instructive to 
briefly discuss an alternative approach J^] to implement- 
ing the gaussian assumption for the field m. Consider a 
given mapping <f>(m) between the original field and the 
auxiliary field m. The spatial variation of <f> across an in- 
terface is asymptotically identical (in the scaling regime) 
to that of an equilibrium interface. It follows that, at 
an interface, |V<^| is a given constant, determined by the 
equilibrium interface profile, and therefore |Vm| is a con- 
stant, which we shall call a. 

This discussion shows that the gaussian approximation 
cannot be exact because, for a gaussian field m, |Vm| 
fluctuates with position on the interface. An exception 
is the limit d — > 00, when |Vm| approaches the limit 
[J2i=i((^ m / 'dxi) 2 )} 1 / 2 . We shall return to this limit be- 
low. 

Replacing |Vm| by a in (||) and (|ll"| ) leads to the prob- 
ability distribution 



P v (v) 



(5 (v- 



S(m)) 



(6(m)) 



Sa 2 \ 1/2 
2nD J 



exp 



Sa 2 
2D 



(17) 



This new approach, therefore, fails to reproduce the 
power-law tail in P v (v) predicted using the general argu- 
ments presented in section II, but gives instead a gaussian 
distribution for all d. While both approaches would seem 
to be equally valid a priori, the calculation of P v (v) dis- 
criminates clearly in favor of the first approach of treating 
the field m consistently as gaussian throughout. In the 
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limit d — > oo, however, both approaches agree, since jl5| ) 
has, for large d, the limiting gaussian form 



P v (v) = (i/ir) 1 ' 2 exp(-tv 2 ) 



(18) 



using v — » (d/2t) 1 ' 2 from the OJK theory. 

A concrete realization of the new scheme requires a 
model which ensures that the interface thickness is time- 
independent. Although it is difficult to enforce the con- 
straint |Vm| — a at interfaces, models in which the 



looser condition ((Vm) 2 



a is imposed can be de- 



vised. The simplest such model is that proposed by Bray 
and Humayun IHJlOfl, equivalent to the Oono-Puri exten- 
sion |ll| of the OJK model. Within this model, the field 
is still strictly gaussian, so P v (v) still has the form (15). 
This model allows, however, a quantitative comparison 
of (0) and © for d -> oo. One finds §,|l§ S = At/d, 
T = (d + 2)/2dt, and U = 1/d (corresponding to a = 1), 
and hence D = 2/d. Putting these results into either 
(15), with d — > oo, or into fll7|), with a = 1, reproduces 

S 



IV. VECTOR FIELDS 

The scaling arguments of section II can be general- 
ized rather simply to n-component vector fields in spa- 
tial dimension d > n (the latter condition being neces- 
sary for the existence of localized defects). The quantita- 
tive (though approximate) methods of section III are, for 
technical reasons, less straightforward to extend to the 
general case. Mazenko Qj has recently given the result 
for nonconserved fields with point defects (n = d). 

The scaling argument can be presented in a rather gen- 
eral way. If the dynamical exponent is z, i.e. L(t) ~ 
i 1 / 2 , it follows that the characteristic velocity scales as 
dL/dt ~ L~( z ~ x \ The next step is to assume that a 
small defect structure [a small domain (n — 1), a small 
vortex loop (n = 2, d = 3) or vortex-antivortex pair 
(n = 2 = d), etc.], of size much less than L(t), collapses 
under its own internal forces at a rate consistent with 
scaling. That is, the scale r of the structure evolves as 
dr/dt ~ — r -( z - 1 ) for r <C L(t). The collapse time of 
such a structure will then scale as r 2 , and structures with 
r < (At) 1 / 2 will disappear in a time interval At. This 
scenario agrees precisely with what we found explicitly 
for scalar fields in section II, where z = 2 (3) for non- 
conserved (conserved) fields. The case n = 2 has been 
discussed by Rutenberg and Bray Jig ]. 

The number of defect structures per unit volume with 
size between r and r + dr is again given by the scaling 
form (|l|) . Requiring that the number of structures which 
disappear in interval At be linear in At forces the scal- 
ing function f{x) in (fil) have the small- argument form 
f(x) ~ x z ~ x , giving n(r) ~ r 2-1 /L d+Z for r -C L. The 
core volume of a structure of size r scales as r d ~™, while 
the total core volume per unit volume of space scales as 



L n . The core- volume weighted probability for defect 
structure size between r and r + dr is, therefore, 



z—l+d—n 



P r [f) dr ~ r 

a generalization of (|] 
to the power-law tail 

Pv{v) - 



dr/L 



d-\-z — n 



r < L 



(19) 



and 



1 



Using v ~ r ( z ^ leads 



(20) 



with tail-exponent 

p = 2 + (d+l-n)/(z- 1) . 



(21) 



This result generalizes (||) and (||), to which it reduces 
for n = 1 and z — 2, 3 respectively. 

We can compare this result with the analytical result 
for P v (v) obtained using the gaussian closure approxi- 
mation for nonconserved fields (z = 2) with point de- 
fects (n = d) M. For this case ( pl| ) becomes p — 3, i.e. 
P v (v) ~ l/L 2 ir. Note that P v (v) is the probability dis- 
tribution for the magnitude of the velocity. The gaussian 
closure calculation gives an expression for the probability 
distribution P(v) of the vector velocity, in the form H 



P(v) 



i 



(7TU 2 )"/ 2 (1 + i; 2 /v 2 )(™+ 2 )/ 2 



(22) 



The normalization of this distribution is J d n v P(v) = 1, 
whereas P v (v) is normalized as dv P v (v) = 1. The re- 
lation between these two distributions is simply P v (v) = 
[2n n / 2 /T(n/2)]v n - 1 P(v), i.e. 



Pv(v) 



(v/v) r 



v (l+-y 2 /« 2 ) (n+2)/2 



(23) 



The power-law tail is P v (v) ~ nv 2 /v 3 , which agrees with 
the scaling result (El]) (for the case n = d, z = 2). 



V. CONCLUSION 

The velocity distribution of topological defects in a 
phase-ordering system has been discussed. A simple scal- 
ing argument, introduced in the context of domain walls 
in section II, and generalized to vector fields in section 
IV, predicts a power-law tail at large velocity. The tail 
exponent p is given by the general result (pl|), where z 
is the dynamical exponent that describes the coarsening 
dynamics, via L(t) ~ t x l z . 

For the special case of non-conserved fields, approxi- 
mate analytical calculations, based on a gaussian closure 
scheme, have been performed for scalar fields (section 
III) and for vector fields with point defects j|. In both 
cases the approximate result exhibits the power-law tail 
predicted by the scaling arguments. Comparison of equa- 
tions (15) and ( |23| ) suggests the following conjecture for 
general nonconserved vector fields, within the gaussian 
approximation: 



4 



Pv(v) 



2r( rf±2 ) 



r(5)r(i 



-n+2 



1 (v/v) 11 - 1 

) i (l + W 2 /w 2 )( d + 2 )/ 2 



(24) 



This result reduces to (|5|) for n = 1 [Lj, to fl23|) for 
?i = d, and gives rise to the power-law tail P v {v) ~ 
with p = rf + 3 — n, in agreement with the scaling predic- 
tion @ for 2 = 2. 

As a final comment we note that for certain systems 
the growth law for L{t) is expected to contain logarith- 
mic corrections to a simple power law. Two examples are 
the nonconserved n = 2 model for d — 2, where one ex- 
pects L{t) ~ (i/lnt) 1 / 2 , and the conserved n = 2 
model for d > 3, where [G^pl ~ (tint) 1 / 4 . A sim- 



ple scaling approach of the type used in sections II and 
IV cannot be applied naively, due to the appearance of a 
new length scale - the defect core size - which enters as 
a short-scale cut-off in the logarithms. We would expect, 
however, that any resulting modifications of the final re- 
sult would be limited to possible logarithmic corrections 
to the power-law tail in P v (v). The dominant power-law 
part will be obtained by inserting the appropriate value 
of z (2 or 4 respectively in the cases discussed above) in 
©■ 
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